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Nanodiamond dust and the energy distribution of quasars 



L. Binette 1 , A. C. Andersen 2 , H. Mutschke 3 , and S. Haro-Corzo 1 



in 
o 
o 

(N 

a* 

(N 



> 

in 

m 

On 
O 

in 
o 

Or 

6 
to 

CO 

c3 



13 



Instituto de Astronomfa, Universidad Nacional Autonoma de Mexico, Apartado Postal 70-264, 04510 Mexico, DF, Mexico 

2 Dark Cosmology Center, Juliane Maries Vej 30, DK-2100 Copenhagen, Denmark 

3 Astrophysikalisches Institut und Universitats-Sternwarte (AIU), SchillergaBchen 3, D-07745 Jena, Germany 



the date of receipt and acceptance should be inserted later 



Abstract. The spectral energy distribution of quasars shows a sharp steepening of the continuum shortward of ~ 1100 A. The 
steepening could be a result of dust absorption. We present a dust extinction model which considers crystalline carbon grains 
and compare it with SMC-like dust extinction consisting of a mixture of silicate grains with graphite or amorphous carbon 
grains. We show that the sharp break seen in individual quasar spectra of intermediate redshif ~ 1-2 can be reproduced 
by dust absorption provided the extinction curve consists of nanodiamonds, composed of terrestial cubic diamonds or of 
diamonds similar to the presolar nanodiamonds found in primitive meteorites. 
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1. Introduction 



The ultraviolet energy distribution of quasars is characterized 
by the so-called "big blue bump", which peaks in vF v at ap- 
proximately 1000 A. The quasar 'composite' spectral energy 
distribution (SED) of Telfer etal. (2002, hereafter TZ02) ob- 
tained by co-adding 332 HST-FOS archived spectra of 184 
quasars between redshifts 0.33 and 3.6, exhibits a steepen- 
ing of the continuum at ~ 1100 A. A fit of this compos- 
ite SED using a broken powerlaw reveals that the powerlaw 
index changes from —0.69 in the near-UV to —1.76 in the 
far-UV (see TZ02). We label this observed sharp steepen- 
ing the 'far-UV break' . In these proceedings, we compare the 
dust absorption that results from nanodiamonds with that of 
more traditional dust species. In a previous proceedings (Bi- 
nette etal. 2005b), we described how we came to consider 
the possibility of carbon crystallite dust. The argumentation 
in support of the dust absorption interpretation of the UV 
break has been fully described in Binette etal. (2005a, here- 
after BM05). Complementary information about crystalline 
dust can be found in this and other proceedings (Binette et al. 
2005b, 2005c). 
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2. Dust extinction model 

Nanodiamonds are the most abundant presolar grain type 
found in the relatively unprocessed meteorites, called car- 
bonaceous chondrites. The diamonds account for about 3% 
of the total amount of carbon in this type of meteorites! The 
main difference between meteoritic diamonds and cubic ter- 
restrial diamonds is that meteoritic diamonds possess surface 
impurities that significantly alter their optical properties. This 
can be appreciated in Fig.^where we compare the extinction 
cross-section resulting from cubic terrestrial diamonds (black 
continuous line) with that from nanodiamonds from the Al- 
lende meteorite (gray continuous line). In both dust models, 
we assume a powerlaw grain size distribution of index —3.5 
bounded by size limits of 2 and 25 A. We assume the grains to 
be spherical and adopt the Mie theory (Bohren and Huffman 
1983) and the complex refraction indices n + ik of Mutschke 
etal. (2004) (for the Allende nanodiamonds) and Edwards & 
Philipp (1985) (for the terrestrial diamonds). It turns out that 
within the above small grain size regime, using a log-normal 
distribution or altering the above grain size limits slightly, do 
not produce any differences in the derived extinction curves, 
as explained in BM05. 

It is apparent from Fig.^that the steepness of the near- 
UV extinction rise in the case of nanodiamonds ought to pro- 
duce a sharp absorption break, a desired feature for any dust 
model that aims at fitting the UV quasar break. The nanodi- 
amond extinction curves in Fig. ^ were normalized in such a 
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Fig. 1. Extinction cross section for terrestrial cubic nanodia- 
monds (black continuous line labeled Dl) and meteoritic nan- 
odiamonds from the Allende Meteorite (gray continuous line 
labeled Al). The silver continuous line represents a graphite 
and silicate model of the ISM dust by Martin & Rouleau 
(1991), with grain sizes ranging from 50 to 2500 A. The dot- 
ted line corresponds to a similar (ISM) grain composition 
but with sizes confined to the range 50 to 300 A. Notice the 
graphite absorption feature at 2175 A in both currves. The 
long dashed line describes a dust model consisting of grains 
composed of silicate (50%), amorphous carbon (45%) and 
5% of graphite, with grain sizes within the range 50 to 300 A. 
The dot-dashed gray line represents the SMC dust model of 
Pei (1992) multiplied here by the factor 5.5. 

way that they represent the case of having all the carbon in 
the dust (assuming a solar C/H abundance ratio). This nor- 
malization would need to be scaled according to the actual 
but unknown dust-to-gas ratio appropriate to the quasar envi- 
ronment. 

3. Powerlaws absorbed by nanodiamond dust 

In BM05, we present calculations in which the dust is either 
intrinsic to the quasars or intergalactic. In the end, we con- 
cluded that only the intrinsic dust hypothesis was satisfactory. 
The calculations presented here will assume the intrinsic dust 
case, which requires gas columns of order 10 20 cm~ 2 assum- 
ing Solar carbon abundance and full depletion onto nanodia- 
mond grains. 

Both nanodiamond types can produce a sharp UV break 
as shown in Fig.|2]where we assumed a powerlaw energy dis- 
tribution of spectral index —0.40 (dotted line) and a column 
of 10 20 cm~ 2 . The opacity is given by r^ xt = N-na^ xt {\) 
and the transmission function is simply e~ T * It is inter- 
esting to notice that the break in the absorbed powerlaws in 
Fig. [2] is shortward (longward) of Lya in the case of cubic 
diamonds (meteoritic nanodiamonds), respectively. 
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Fig. 2. Absorption models assuming a powerlaw SED of in- 
dex a u — —0.4 (dotted line). The nanodiamond dust screen 
corresponds to a hydrogen column of 10 20 cm~ 2 . The ab- 
sorbed SED represented by the continuous line corresponds 
to the extinction by terrestrial cubic diamonds (Dl), while 
the long dashed line corresponds to nanodiamonds from the 
Allende meteorite (Al). The gray lines correspond to a mix- 
ture of the two flavors: dashed line: 60% Dl + 40% Al, dot- 
dashed line: 30% Dl + 70% Al. Labelled pointers indicate 
where emission lines are expected in the rest-frame quasar 
spectra. 

It turns out that the UV break position in most spectra 
is about midway between the break positions produced by 
either of the two nanodiamond types, suggesting that a com- 
bination of both types is needed, as described in more detail 
in BM05. We found that most quasar spectra favor a com- 
position consisting on the order of 30% of cubic diamonds 
and 70% of meteoritic nanodiamonds. In Fig.[3k, we illus- 
trate the effect of varying the absorption column A r n, assum- 
ing the above mixture (30% and 70%, respectively, equiva- 
lent to fni = 0.3). The curves shown correspond to the four 
H columns values of ./V20 = 0.5, 1.0, 2.0 and 4.0, in units of 
10 20 cm~ 2 . The same powerlaw SED (dotted line) is assumed 
as for Fig.|2] that is ct v = —0.4. 

4. Comparison with SMC-like dust 

Can ISM or SMC-like dust reproduce the quasar break? Ac- 
cording to Shang et al. (2005), reddening by ISM or SMC- 
like grains "is not able to produce the spectral break seen in 
the AGN sample, without leaving a clear signature at longer 
wavelengths". Instead of a sharp break, ISM extinction re- 
sults in a shallow rollover that extends in the near-UV and 
even in the optical domain. In the HST-FOS sample that we 
studied, most spectra (the 50 class A and B spectra) are quite 
hard in the near-UV, consistent with a powerlaw continuum 
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Fig. 3. Absorbed powerlaw SED of index a v = —0.4 (dotted line) for the following Ah columns in units of 10 20 cm -2 : 
A20 = 0.5 (dashed line), 1.0 (solid line), 2.0 (long dashed line) and 4.0 (dash-dot line). The following two dust compositions 
are considered. Panel a (left): the dust corresponds to a combination of terrestrial cubic diamonds (30%) and primitive mete- 
orite nanodiamonds (70%). Panel b (right): the dust grains are made of silicate (50%), amorphous carbon (45%) and graphite 
(5%). 
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Fig. 4. The rest-frame spectrum of class (A) quasar 
PG 1008+1319. Notice the far-UV rise shortward of the UV 
break. The nanodiamond dust extinction model (gray con- 
tinuous line) assumes a powerlaw with a v = +0.13 (from 
Neugebauer etal. 1987), A20 = 1.2 and the same dust com- 
position as in Fig.|5t (mixed nanodiamond types). The gray 
dashed line corresponds to extinction by silicate-amorphous 
carbon dust (as in Fig.[3p) assuming a column A20 = 2.4. 
The transmitted flux in the latter model has been multiplied 
by 1.3 so that longward of 1800 A it superimposes the trans- 
mitted flux of the nanodiamond model. 



longward of 1 100 A. They do not indicate any substantial red- 
dening longward of the break (except the 8 objects associated 
to class C, see BM05). In Fig.^ the extinction cross-section 
resulting from ISM dust is represented by the silver continu- 
ous line. The model depicted consists of silicate and graphite 
grains with grain sizes ranging between 50 and 2500 A. 

The extinction by SMC-like dust, however, is more peaky 
in the near-UV and, as will be shown, cannot reproduce the 
sharp 1000 A break either. For comparison, the long dashed 
line in Fig. represents the cross-section of a dust model 
consisting of relatively small grain sizes, ranging from 3 to 
300 A, but with graphite replaced by amorphous carbon (us- 
ing results from Rouleau & Martin 1991), that is, dust com- 
posed of 50% silicate grains, 45% amorphous C grains and 
5% graphite grains. The dust-to-gas ratio is the same as for 
the Galactic ISM model. The reason for replacing most of 
the 'usual' graphite grains by amorphous carbon is guided 
by the need to remove the 2175 A absorption feature present 
when small graphite grains are abundant. Our aim is to test 
an extinction curve that would be SMC-like, in accordance 
with the extinction law inferred from the SLOAN quasar 
sample by Richards etal. (2003), Hopkins etal. (2004) and 
Willott (2005). In Fig|3] we indicate how the previous power- 
law SED would look if absorbed by dust containing silicates 
and amorphous carbon (corresponding to the long dashed 
line extinction model in Fig. 1), for the four columns val- 
ues of A 20 = 0.5, 1.0, 2.0 and 4.0, in units of 10 20 cm- 2 . 
A sharp break does occur for columns A20 ~ 2. The break 
takes place, however, at too long a wavelength in the near-UV. 
This is illustrated in Fig|4] where we compare the spectra of 
PG 1008+1319 with a powerlaw absorbed by nanodiamonds 
(gray continuous line) and by silicate-amorphous carbon dust 
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(gray dash line). Clearly, only the nanodiamond dust provides 
a satisfactory fit. 

Adopting the more commonly used SMC dust model of 
Pei (1992), which consists only of silicate grains, would not 
improve the situation, since the extinction cross-section is 
even shallower in the near-UV in the Pei model than in the 
case of the silicate-amorphous carbon dust. A comparison 
of these two extinction curves is made in Fig.^ The Pei 
cross-section has been mutiplied by 5.5 Fig.[0 to show that 
it approximately overlaps in the range 1800-5000 A with the 
silicate-amorphous C extinction curve. This confirms that the 
selective extinction of both dust models is quite comparable 
in the optical domain. 

5. Conclusion 

Dust absorption by nanodiamonds is successful in reproduc- 
ing the 1000 A break as well as the far-UV rise seen at shorter 
wavelengths in distant quasars. To confirm the possible exis- 
tence of nanodiamond grains in quasars, observations of the 
far-infrared emission bands caused by hydrogenated nanodi- 
amonds (van Kerckhoven et al. 2002; Jones et al. 2004) could 
be attempted or, alternatively, one may try to confirm a flux 
rise shortward of 700 A for as many quasars as possible. The 
latter would require high quality far-UV sensitive observa- 
tions, which will hopefully be provided by the World Space 
Observatory satellite (Barstow etal. 2003). 

The model presented is successful in both fitting the 
breaks shape and the position of, as well as in reproducing the 
sharp flux recovery observed in key quasar spectra, around 
660 A. While very high redshift quasars (z J; 2.5) appear not 
to be absorbed by dust, their energy distribution suggests the 
existence of a higher energy continuum rollover at ~ 670 A. 
As shown in BM05, when such a break is incorporated into 
the description of the energy distribution of all quasars, the 
dust absorption model can account for the overall shape of the 
'composite' spectral energy distribution of the 184 quasars 
constructed by TZ02. 
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